Isostaticity, auxetic response, surface modes, and conformal invariance in twisted 
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Model lattices consisting of balls connected by central-force springs provide much of our under- 
standing of mechanical response and phonon structure of real materials. Their stability depends 
critically on their coordination number z. d-dimensional lattices with z = 2d are at the threshold of 
mechanical stability and are isostatic. Lattices with z < 2d exhibit zero-frequency "floppy" modes 
that provide avenues for lattice collapse. The physics of systems as diverse as architectural struc- 
tures, network glasses, randomly packed spheres, and biopolymer networks is strongly influenced by 
a nearby isostatic lattice. We explore elasticity and phonons of a special class of two-dimensional 
isostatic lattices constructed by distorting the kagome lattice. We show that the phonon structure of 
these lattices, characterized by vanishing bulk moduli and thus negative Poisson ratios and auxetic 
elasticity, depends sensitively on boundary conditions and on the nature of the kagome distortions. 
We construct lattices that under free boundary conditions exhibit surface floppy modes only or a 
combination of both surface and bulk floppy modes; and we show that bulk floppy modes present 
under free boundary conditions are also present under periodic boundary conditions but that sur- 
face modes are not. In the the long-wavelength limit, the elastic theory of all these lattices is a 
conformally invariant field theory with holographic properties, and the surface waves are Rayleigh 
waves. We discuss our results in relation to recent work on jammed systems. Our results highlight 
the importance of network architecture in determining floppy-mode structure. 

PACS numbers: 



O 

o 



> 
o 



X 



I. INTRODUCTION 

Networks of balls and springs or frames of nodes con- 
nected by compressible struts provide realistic models for 
physical systems from bridges to condensed solids. Their 
elastic properties depend on their coordination number z 
- the average number of nodes each node is connected to. 
If z is large enough, the networks are elastic solids whose 
long- wavelength mechanical properties are described by a 
continuum clastic energy with non- vanishing elastic mod- 
uli. If z is small enough, the networks have deformation 
modes of zero energy - they are floppy. As z is increased 
from the floppy side, a critical value, z c , is reached at 
which springs provide just enough constraints that the 
system has no zero-energy "floppy" modes [l[ (or mech- 
anisms @ in the engineering literature), and the sys- 
tem is isostatic. The phenomenon of rigidity percolation 
[H, 0] whereby a sample spanning rigid cluster develops 
upon the addition of springs is one version of this floppy- 
to-rigid transition. The coordination numbers of whole 
classes of systems, including engineering structures 0, H[ 
(bridges and buildings), randomly packed spheres near 
jamming &ticl : network glasses [l|,[llj], cristobalites [l2j], 
zeolites (l3l Il4|. and biopolymer networks [T5l - [l8j are 
close enough to z c that their elasticity and mode struc- 
ture is strongly influenced by those of the isostatic lattice. 

Though the isostatic point always separates rigid from 
floppy behavior, the properties of isostatic lattices are 
not universal; rather they depend on lattice architecture. 
Here we explore the the unusual properties of a partic- 
ular class of periodic isostatic lattices derived from the 



two-dimensional kagome lattice by rigidly rotating trian- 
gles through an angle a without changing bond lengths 
as shown in Fig.[T] The bulk modulus B of these lattices 
is rigorously zero for all a ^ 0. As a result, their Poisson 
ratio acquires its limit value of —1; when stretched in 
one direction, they expand by an equal amount in the or - 
thogonal direction: they are maximally auxetic 
These modes represent collapse pathways [H, [24 
kagome lattice. Modes of isostatic systems are generally 
very sensitive to boundary conditions [2oT - |27| , but the de- 
gree of sensitivity depends on the details of lattice struc- 
ture. For reasons we will discuss more fully below, modes 
of the square lattice, which is isostatic, are in fact insen- 
sitive to changes from free boundary conditions (FBCs) 
to periodic boundary conditions (PBCs), whereas those 
of the undistorted kagome lattice are only mildly so. The 
modes of both, however, change significantly when rigid 
boundary conditions (RGBs) are applied. We show here 
that in all families of the twisted kagome lattice, modes 
depend sensitively on whether FBCs, PCBs or RGBs are 
applied: finite lattices with free boundaries have floppy 
surface modes that are not present in their periodic or 
rigid spectrum or in that of finite undistorted kagome 
lattices. In the long wavelength limit, the surface floppy 
modes, which are present in any 2d material with B = 0, 
reduce to surface Rayleigh waves [28| described by a con- 
formally invariant energy whose analytic eigenfunctions 
are fully determined by boundary conditions. At shorter 
wavelengths, the surface waves become sensitive to lattice 
structure and remain confined to within a distance of the 
surface that diverges as the undistorted kagome lattice is 
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FIG. 1: (a) Section of a kagome lattice with N x = N y = 4 and 
N c — N x N y three-site unit cells. Nearest-neighbor bonds, oc- 
cupied by harmonic springs, are of length a. The rotated row 
(second row from the top) represents a floppy mode. Next- 
nearest neighbor bonds are shown as dotted lines in the lower 
left hexagon. The vectors ei, e2, and e3 indicate symmetry 
directions of the lattice. The numbers in the triangles indicate 
those that twist together under PBCs in zero modes along the 
three symmetry direction. Note that there are only 4 of these 
modes, (b) Section of a square lattice depicting a floppy mode 
in which all sites along a line are displaced uniformly, (c) 
Twisted kagome lattice, with lattice constant at = 2a cos a, 
derived from the undistorted lattice by rigidly rotating tri- 
angles through an angle a. A unit cell, bounded by dashed 
lines, is shown in violet. Arrows depict site displacements 
for the zone-center, i.e., zero wavenumber, <f> mode which has 
zero (nonzero) frequency under free (periodic) boundary con- 
ditions. Sites 1, 2, and 3 undergo no collective rotation about 
their center of mass whereas sites 1, 2*, and 3* do. (d) Su- 
perposed snapshots of the twisted lattice showing decreasing 
areas with increasing a. 



approached. In the simplest twisted kagome lattice, all 
floppy modes are surface modes, but in more complicated 
lattices, including ones with uniaxial symmetry, we con- 
struct, there are both surface and bulk floppy modes. 

Arguments due to J.C. Maxwell [29] provide a criterion 
for network stability: networks in d dimensions consist- 
ing of TV nodes, each connected with central-force springs 
to an average of z neighbors, have Nq = dN — \zN zero- 
energy modes when z < 2d (in the absence of redundant 
bonds - see below). Of these a number, TV tr , which de- 



pends on boundary conditions, are trivial rigid transla- 
tions and rotations, and the and the remainder are floppy 
modes of internal structural rearrangement. Under FBCs 
an PBCs, TV tr equals d(d+l)/2 and d, respectively. With 
increasing z, mechanical stability is reached at the iso- 
static point at which Nq = TV tr . The Maxwell argument 
is a global one; it does not provide information about 
the nature of the floppy modes and does not distinguish 
between bulk or surface modes. 



II. KAGOME ZERO MODES AND ELASTICITY 

The kagome lattice of central force springs shown in 
Fig. [Tfa) is one of many locally isostatic lattices, includ- 
ing the familiar square lattice lattice in two dimensions 
[Fig. QJb)] and the cubic and pyrochlore lattices in three 
dimensions, with exactly z = 2d nearest-neighbor (TV TV) 
bonds connected to each site not at a boundary. Un- 
der PBCs, there are no boundaries, and every site has 
exactly 2d neighbors. Finite, TV-site sections of these 
lattices have surface sites with fewer than 2d neighbors 
and of order vTV zero modes. The free kagome lattice 
with N x and N y unit cells along its sides [Fig. [Ha)] nas 
TV = 3TV x TVj, sites, TV B = 6N x N y - 2(N X + N y ) + 1 bonds, 
and TVo = 2(N X + TV y ) — 1 zero modes, all but three 
of which are floppy modes. These modes, depicted in 
Fig-Ufa), consist of coordinated counter rotations of pairs 
of triangles along the symmetry axes ei , e2 and e3 of the 
lattice. There are N x modes associated with lines par- 
allel to ei, N y associated with lines parallel to e3, and 
TVr + N y — 1 modes associated with lines parallel to e2. 

In spite of the large number of floppy modes in the 
kagome lattice, its longitudinal and shear Lame coeffi- 
cients, A and /i, and its Bulk modulus B = A + /i are 
nonzero and proportional to the nearest neighbor (TV TV) 
spring constant k: A = [i = ^/3k/8, and B = A + fi = 
-\/3fc/4. The zero modes of this lattice can be used to 
generate an infinite number of distorted lattices with un- 
stretched springs and thus zero energy [H], [3(| ■ We con- 
sider only periodic lattices, the simplest of which are the 
twisted kagome lattices obtained by rotating triangles of 
the kagome unit cell through an angle a as shown in 
Figs.Hfc) and (d) (MIHI- These lattices have C 3v rather 
than Cq v symmetry and, like the undistorted kagome lat- 
tice, three sites per unit cell. As Fig. [T^d) shows, the 
lattice constant of these lattices is cll = 2acosev, and 
their area A a decreases as cos 2 a as a increases. The 
maximum value that a can achieve without bond cross- 
ings is 7r/3 so that the maximum relative area change 
is A^/s/Aq = 1/4. Since all springs maintain their rest 
length, there is no energy cost for changing a, and as a 
result, B is zero for every a ^ 0, whereas the shear mod- 
ulus fj, — V3k/8 remains nonzero and unchanged. Thus, 
the Poisson ratio a = (B — fi)/(B + /i) attains its small- 
est possible value of —1. For any a ^ 0, the addition of 
next-nearest-neighbor (TV TV TV) springs, with spring con- 
stant k' (or of bending forces between springs) stabilizes 
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FIG. 2: phase diagram in the a — k' plane showing region 
with negative Poisson ratio a. 



zero-frequency modes and increases B and a. Neverthe- 
less, for sufficiently small k', a remains negative. Figure 
[2] shows the region in the k! — a plane with negative a. 
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FIG. 3: (a)Phonon spectrum for the undistorted kagome lat- 
tice. Dashed lines depict frequencies at k — and full lines 
at k' > 0. The inset shows the Brillouin zone with symme- 
try points F, M, and K. Note the line of zero modes along 
YM when k! = 0, all of which develop nonzero frequencies 
for wavenumber q > when k' > reaching oj* ~ y/W on 
a plateau beginning at q ~ q* ~ \fk! defining a length scale 
V — 1/q*. (b) Phonon spectrum for a > and k' — 0. The 
plateau along YM defines uj a ~ vfe|sin.a| and its onset at 
q a ~ Lu a defines a length l a ~ f/| sina|. 



III. KAGOME PHONON SPECTRUM 

We now turn to the linearized phonon spectrum of the 
kagome and twisted kagome lattices subjected to PBCs. 
These conditions require displacements at opposite ends 
of the sample to be identical and thus prohibit distortions 
of the shape and size of the unit cell and rotations but 
not uniform translations, leaving two rather than three 
trivial zero modes. The spectrum [32[ of the three lowest 
frequency modes along symmetry directions of the undis- 
torted kagome lattice with and without NNN springs 
is shown in Fig. [U[a). When k' — 0, there is a floppy 
mode for each wavenumber q ^ running along the 
entire length of the three symmetry-equivalent straight 
lines running from M to T to M in the Brillouin zone 
[See inset to Fig. [3]- When N x = N y , there are exactly 
N x — 1 wavenumbers with q ^ along each of these 
lines for a total of 3(N X — 1) floppy modes. In addition, 
there are three zero modes at q = corresponding to 
two rigid translations, and one floppy mode that changes 
unit cell area at second but not first order in displace- 
ments, yielding a total of 3N X zero modes rather than 
the AN X — 1 modes expected from the Maxwell count un- 
der FBCs. This is our first indication of the importance 
of boundary conditions. The addition of NNN springs 
endows the floppy modes at k' = with a characteris- 
tic frequency w* ~ \fk! and causes them to hybridize 
with the acoustic phonon modes [Fig. EJa)] [32]. The 
result is an isotropic phonon spectrum up to wavenum- 
ber q* — 1/7* ~ -v/fc 7 and gaps at T and M of order w*. 
Remarkably, at nonzero a and k' — 0, the mode struc- 
ture is almost identical to that at a = and k' > 
with characteristic frequency ui a ~ y/k\ sina| and length 
l a ~ l/uj a . In other words, twisting the kagome lattice 
through an angle a has essentially the same effect on the 
spectrum as adding NNN springs with spring constant 
| sina| 2 fc. Thus under PBCs, the twisted kagome lattice 
has no zero modes other than the trivial ones: it is "col- 



lectively" jammed in the language of references (27l. [33|. 
but because it is not rigid with respect to changing the 
unit cell size, it is not strictly jammed. 



IV. MODE COUNTING AND STATES OF SELF 
STRESS 

To understand the origin of the differences in the zero- 
mode count for different boundary conditions, we turn to 
an elegant formulation Q of the Maxwell rule that takes 
into account the existence of redundant bonds (i.e., bonds 
whose removal does not increase the number of floppy 
modes and states in which springs can be under states 
of self-stress. Consider a ring network in two dimensions 
shown in Fig. 2] with N = 4 nodes and 2V& = 4 springs 
with three springs of length a and one spring of length 
b. The Maxwell count yields ^0 = 4 = 3 + 1 zero modes: 
two rigid translations, one rigid rotation, and one internal 
floppy mode - all of which are "finite-amplitude" modes 
with zero energy even for finite-amplitude displacements. 
When b = 3a, the Maxwell rule breaks down. In the zero- 
energy configuration, the long spring and the three short 
ones are colinear, and a prestressed state in which the b- 
spring is under compression and the three a-springs are 
under tension (or vice versa) but the total force on each 
node remains zero becomes possible. This is called a state 
of self-stress. The system still has three finite amplitude 
zero modes corresponding to arbitrary rigid translations 
and rotations, but the finite-amplitude floppy mode has 
disappeared. In the absence of prestress, it is replaced 
by two "infinitesimal" floppy modes of displacements of 
the two internal nodes perpendicular of the now linear 
network. In the presence of prestress, these two modes 
have a frequency proportional to the square root of the 
tension in the springs. Thus, the system now has one 
state of self stress and one extra zero mode in the absence 



of prestress, implying No = 2N — Nb + S, where S is the 
number of states of self stress. 

This simple count is more generally valid as can be 
shown with the aid of the equilibrium and compatibility 
matrices @, denoted, respectively, as H and C = H T . H 
relates the vector t of Nb spring tensions to the vector 
f of dN forces at nodes via H ■ t = f , and C relates 
the the vector d of dN node displacements to the vector 
e of Nb spring stretches via C ■ d = e. The dynamical 
matrix determining the phonon spectrum is D = fcHH T . 
Vectors to in the null space of H, (H - to =0), describe 
states of self-stress whereas vectors do in the null space of 
C represent displacements with no stretch e, i.e., modes 
of zero energy. Thus the nullspace dimensions of H and 
C are, respectively, S and Nq. The rank- nullity theorem 
of linear algebra 34j states that the rank r of a matrix 
plus the dimension of its null space equals its column 
number. Since the rank of a matrix and its transpose 
are equal, the H and C matrices, respectively, yield the 
relations r + S = Nb and r + Nq — dN, implying Nq = 
dN — Nb + S. Under PBCs, locally isostatic lattices 
have z = 2d exactly, and the Maxwell rule yields No = 
0: there should be no zero modes at all. But we have 
just seen that both the square and undistorted kagome 
lattices under PBCs have of order y/N zero modes as 
calculated from the dynamical matrix, which, because it 
is derived from a harmonic theory, does not distinguish 
between infinitesimal and finite-amplitude zero modes. 
Thus, in order for there to be zero modes, there must be 
states of self-stress, in fact one state of self-stress for each 
zero mode. 

In the square lattice under FBCs, N = N x N y and 
Nb = 2N x N y — N x — N y , there are no states of self stress, 
and Nq = N x + N y zero modes depicted in Fig. QJb). 
Under PBCs, the dimension of the nullspace of H is 
S = N x + N y , and there are also N = S = N x + N y 
zero modes that are identical to those under FBCs. We 
have already seen that there are No = 2(N X + N v ) — 1 
zero modes in the free undistorted kagome lattice. Di- 
rect evaluations [23| (See Text SI) of the dimension of 
the null spaces of H and C for the undistorted kagome 
lattice with PBCs yields S = N = 3N X when N x = N y . 
The zero modes under PBCs are identical to those un- 
der FBCs except that the 2N X — 1 modes associated with 
lines parallel to e2 under FBCs get reduced to N x modes 
because of the identification of apposite sides of the lat- 
tice required by the PBCs as shown in Fig. QJa). Thus 
the modes of both the square and kagome lattices do 
not depend strongly on whether FCBs or PBCs are ap- 
plied. Under RBC's, however, the floppy modes of both 
disappear. The situation for the twisted kagome lattice 
is different. There are still 2(N X + N y ) — 1 zero modes 
under FBCs, but there are only two states of self stress 
under PBCs and thus only N = S = 2 zero modes, as 
a direct evaluation of the null spaces of H and C verifies 
(See Appendix for details) , in agreement with the results 
obtained via direct evaluation of the eigenvalues of the 
dynamical matrix [32l |35| . All of the floppy modes under 
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FIG. 4: (a) Ring-network with b > 3a showing internal floppy 
mode, (b) Ring-network with b — 3a showing one of the two 
infinitesimal modes. 



FBCs have disappeared. 



V. EFFECTIVE THEORY AND EDGE MODES 

An effective long-wavelength energy E c g for the low- 
energy acoustic phonons and nearly floppy distortions 
provides insight into the nature of the modes of the 
twisted kagome lattice. The variables in this theory are 
the vector displacement field u(x) of nodes at undistorted 
positions x and the scalar field </>(x) describing nearly 
floppy distortions within a unit cell. The detailed form 
of E e ff depends on which three lattice sites are assigned 
to a unit cell. Figure [T](c) depicts the lattice distortion 
<f> for the nearly floppy mode at T (with energy propor- 
tional to | sina| 2 ) along with a particular representation 
of a unit cell, consisting of a central asymmetric hexagon 
and two equilateral triangles, with 8 sites on its bound- 
ary. If sites 1, 2, and 3 are assigned to the unit cell, 
then the distortion <p involves no rotations of these sites 
relative to their center of mass, and the harmonic limit 
of E c g depends only on the symmetrized and linearized 
strain u 



(diUj + djUi)/2 and on <f>: 



E = - 
2 



d 2 x 



2fJLU^ 



+K{cj> + £u u ) 2 + V{di4>f - WY ijk Uij d k (j>\ , (1) 

where = tiy — -^SijUkk is the symmetric-traceless 
stain tensor, fi = y/3k/8, K = 3\/3 tan 2 a/a 2 , £ = 
a esc ql /(2a/3), W = \/3k/4 + 0{a 2 ), and V = \/3k/8 + 
0(a 2 ). The last term in which T t: > k is a third-rank ten- 
sor, whose only non-vanishing components are T xxx = 
_ T xyy = _ V yyx = _ T xyx = is invariant un d er op- 
erations of the group Cq v but not under arbitrary rota- 
tions. The K^4>ua term is the only one that reflects the 
C^ v (rather than Cq v ) symmetry of the lattice. There 
are several comments to make about this energy. The 
gauge-like coupling in which the isotropic strain un ap- 
pears only in the combination (<f> + £uu) 2 guarantees 
that the bulk modulus vanishes: <p wm simply relax to 
—^ua to reduce to zero the energy of any state with 
nonvanishing ua. The coefficient K can be calculated 
directly from the observation that under <fi alone, the 
length of every spring changes by 5a = — \/30sina, and 
this length change is reversed by a homogenous volume 
change ua — SA a /A a — —2Sa/a. In the a —¥ limit, 
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K — > 0, and the energy reduces to that of an isotropic 
solid with bulk modulus B — linic^o K £ 2 = \/3A;/4 if 
the V and W terms, which are higher order in gradients, 
are ignored. The W term gives rise to a term, singular in 
gradients of u, when 4> is integrated out that is respon- 
sible for the deviations of the finite-wavenumber elastic 
energy from isotropy. At small a, the length scale l a ap- 
pears in several places in this energy: in the length £ and 
in the ratios y/fi/K, y/V/K, and y/WjK. At length 
scales much larger than l a , the V and W terms can be 
ignored, and </> relaxes to — £u« leaving only the shear 
elastic energy of an elastic solid proportional fiv^j. At 
length scales shorter than l a , <f> deviates from — £ua and 
contributes significantly to the form of the energy spec- 
trum. If 1, 2* and 3* in Fig.QJd) are assigned to the unit 
cell, then (f> involves rotations relative to the lattice axes, 
and the energy develops a Cosserat-like form [36|,[37j that 
is a function of — a(V x u) z /2 rather than <f>. 

The modes of our elastic energy in the long- wavelength 
limit (ql a <C 1) are simply those of an elastic medium 
with B = 0. In this limit, there are transverse and lon- 
gitudinal bulk sound modes with equal sound velocities 
ct = V ' \ l l ' P = {a/2)y/k/m and cl = \/ (B + fj)/p -> ct, 
where m is the particle mass at each node and p is the 
mass density. In addition there are Rayleigh surface 
waves [28.] in which there is a single decay length (rather 
than the two at B > 0), and displacements are propor- 
tional to e~ qyV cos[q x x] with q y = q x for a semi-infinite 
sample in the right half plane so that the penetration 
depth into the interior is l/q x . These waves have zero 
frequency in two dimensions when B = 0, and they do 
not appear in the spectrum with PBCs. Thus this sim- 
ple continuum limit provides us with an explanation for 
the difference between the spectrum of the free and pe- 
riodic twisted kagome lattices. Under FBCs, there are 
zero-frequency surface modes not present under PBCs. 

Further insight into how boundary conditions affect 
spectrum follows from the observation that the contin- 
uum elastic theory with B = depends only on Uij. 
The metric tensor ffij(x) of the distorted lattice is re- 
lated to the strain Uy (x) via the simple relation gij (x) = 
Sij + 2uy(x); and iiij = [e^(x) - i%g fefe (x)]/2, which is 
zero for g^ = Sij, is invariant, and thus remains equal to 
zero, under conformal transformations that take the met- 
ric tensor from its reference form Sij to h(pt)6ij for any 
continuous function /i(x). The zero modes of the theory 
thus correspond simply to conformal transformations, 
which in two dimensions are best represented by the com- 
plex position and displacement variables z = x + iy and 
w(z) — u x (z) + iuy(z). All conformal transformations 
are described by an analytic displacement field w(z). 
Since by Cauchy's theorem, analytic functions in the in- 
terior of a domain are determined entirely by their val- 
ues on the domain's boundary (the "holographic" prop- 
erty (38|), the zero modes of a given sample are simply 
those analytic functions that satisfy its boundary con- 
ditions. For example, a disc with fixed edges (u = 0) 
has no zero modes because the only analytic function 



satisfying this FBC is the trivial one w(z) — 0; but 
a disc with free edges (stress and thus strain equal to 
zero) has one zero mode for each of the analytic functions 
w(z) — a n z n for integer n > 0. The boundary conditions 
lim :E _ i . 00 u(x, y) = and u(x, y) = u(x + L, y) on a semi- 
infinite cylinder with axis along x are satisfied by the 
function w(z) = e lqxZ = e lqxX e~ qxV when q x = 2nn/L, 
where n is an integer. This solution is identical to that 
for classical Rayleigh waves on the same cylinder. Like 
the Rayleigh theory, the conformal theory puts no restric- 
tion on the value of n (or equivalently q x ). Both theories 
break down, however, at q x — q c ~ min(^ , a -1 ) beyond 
which the full lattice theory, which yields a complex value 
of q y = q' y + iq'y, is needed. 

Figure[5ja) shows an example of a surface wave. At the 
bottom of this figure, u y (x) is an almost perfect sinusoid. 
As y decreases toward the surface, the amplitude grows, 
and in this picture reaches the nonlinear regime by the 
time the surface at y = is reached. Figure OJb) plots 
q' y as a function of q x obtained both by direct numerical 
evaluation and by an analytic transfer matrix procedure 
39] for different values of a (Text SI). The Rayleigh 
limit q' y = q x is reached for all a &s q x ~> 0. Interestingly 
the Rayleigh limit remains a good approximation up to 
values of q x that increase with increasing a. The inset 
to Fig. plots q' y l a as a function of r\ = q x l a and shows 
that in the limit a — > (l a /a — > oo), q' y obeys an a- 
independent scaling law of the form q' — l^ 1 f(q x l a ). The 
full complex q y obeys a similar equation. This type of 
behavior is familiar in critical phenomena where scaling 
occurs when correlation lengths become much larger than 
microscopic lengths. The function f{rf) approaches rj as 
•q — > and asymptotes to 4/3 for 77—5-00. Thus for 
q x l a < 1, qy = q x and for q x l a » 1, q' y = (4/3)/" 1 . 
As a increases, l a /a is no longer much larger than one, 
and deviations from the scaling law result. The situation 
for surfaces along different directions (e.g., along x — 
rather than y = 0) is more complicated and will be 
treated in a future publication [ipj . 



VI. OTHER LATTICES AND RELATION TO 
JAMMING 

The Ca v twisted kagome lattice is the simplest of many 
lattices that can be formed from the kagome and other 
periodic isostatic lattices. Figures [6] (a) and (b) show 
two other examples of isostatic lattices constructed from 
the kagome lattice. Most intriguing is the lattice with 
pgg symmetry. Its geometry has uniaxial symmetry, yet 
its long-wavelength elastic energy is identical to that of 
the Cs v twisted kagome lattice, i.e, it is isotropic with a 
vanishing bulk modulus, and its mode structure near q — 
is isotropic as shown in Fig. [S] (c). Thus, this system 
loses long-wavelength zero-frequency bulk modes of the 
undistorted kagome lattice to surface modes. However, at 
large wavenumber, lattice anisotropy becomes apparent, 
and (infinitesimal) floppy bulk modes appear. Thus in 
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FIG. 5: (a)Lattice distortions for a surface wave on a cylinder, 
showing exponential decay of the surface displacements into 
the bulk. This figure was constructed by specifying a small 
sinusoidal modulation on the bottom boundary and propa- 
gating lattice-site positions upward to the free boundary at 
the top under the constraint of of constant lengths and pe- 
riodic boundary conditions around the cylinder. Distortions 
near and at the top boundary, which have become nonlin- 
ear, are not described by our linearized treatment, (b) q' y CLL 
as a function of q x a,L for lattice Rayleigh surface waves for 
a — 7r/20, 7r/10, 37r/20, n/5, n/A, in order from bottom to top. 
Smooth curves are the analytic results from a transfer matrix 
calculation, and dots are from direct numerical calculations. 
The dashed line is the continuum Rayleigh limit q' y = q x . 
Curves at smaller a break away from this curve at smaller 
values of q y than do those at large a. At a = n/4, q' x di- 
verges at q y a,L = vr. The inset plots q' y l a as a function of 
q x l a for different a. The lower curve in the inset(black) is the 
a- independent scaling function of q y l a reached in the a — s> 
limit. The other curves from top to bottom are for a = n/25, 
7r/12, 7r/9, and n/6 (chosen to best present results rather than 
to match the curves in the main figure). Curves for a < tv/15 
are essentially indistinguishable from the scaling limit. The 
curve at a — n/6 stops because q y < n/aL- 



this and related systems, a fraction of the zero modes of 
the under FBCs are bulk modes that are visible under 
PBCs, and a fraction are surface modes that are not. 

Randomly packed spheres above the jamming tran- 
sition with average coordination number z = 2d + Az 
exhibit a characteristic frequency uj* ~ Az and length 
I* ~ (Az) -1 and a transition from a Debye-like (~ cj d-1 ) 
to a fiat density of states at uj ps uj* [4l|, |42j ■ The square 
and kagome lattices with randomly added NNN springs 
have the sa me p roperties (43l liij . A general "cutting" ar- 
gument [IB), 26] provides a procedure for perturbing away 
from the isostatic limit and an explanation for these prop- 
erties. However, it only applies provided a finite fraction 
of the of order floppy modes of a sample with sides 



FIG. 6: (a) Kagome-based lattice with pgg space group sym- 
metry and uniaxial C2v point group symmetry, (b) Lattice 
with p6 space symmetry but global Ce point-group symme- 
try, (c) Density plot of the spectrum of the lowest frequency 
branch of the pgg uniaxial kagome lattice. The spectrum is 
absolutely isotropic near the origin point F, but it has a zero 
modes on two symmetry related continuous curves at large 
values of wave-number. 

of length L cut from an isostatic lattice with PBCs are 
extended, i.e., have wave functions that extend across the 
sample rather than remaining localized either in the in- 
terior or at the surface the sample. Clearly the twisted 
kagome lattice, whose floppy modes are all surface modes, 
violates this criterion; and indeed, the density of states of 
the lattice with Az = shows Debye-behavior crossing 
over to a flat plateau at uj bs ui a . Adding next nearest 
neighbor bonds gives rise to a length l c (Z" 1 + P -1 ) -1 
and crossover to the plateau at u> c ~ i" 1 . The pgg lat- 
tice in Fig.lHta), however, has both extended and surface 
floppy modes, so its crossover to the a flat plateau occurs 
at uj s» uj* rather than at uj a or uj c . 

VII. CONNECTIONS TO OTHER SYSTEMS 

Our study highlights the rich and remarkable variety 
of physical properties that isostatic systems can exhibit. 
Under FBCs, floppy modes can adopt a variety of forms, 
from all being extended to all being localized near sur- 
faces to a mixture of the two. Under PBCs, the pres- 
ence of floppy modes depends on whether the lattice can 
or cannot support states of self stress. When a lattice 
exhibits a large number of zero-energy edge modes, its 
mechanical/dynamical properties become extremely sen- 
sitive to boundary conditions, much as do the electronic 
properties of the topological states of matter studied in 
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quantum systems (45l - l48j . The zero-energy edge modes 
observed in our isostatic lattices are collective modes 
whose amplitudes decay exponentially from the edge with 
a finite decay length, in direct contrast to the very lo- 
calized and trivial floppy modes arising from dangling 
bonds. We focussed primarily on high-symmetry lat- 
tices derived from the kagome lattice, but the properties 
they exhibit, namely a deficit of floppy modes in the bulk 
and the existence of floppy surface modes, are shared by 
any two-dimensional system with a vanishing bulk mod- 
ulus (or the equivalent in anisotropic systems). Three- 
dimensional analogs of the twisted kagome lattice [3(| can 
be constructed by rotating tetrahedra in pyrochlore and 
zeolite lattices [HI, [3] and in cristobalites [13]. These lat- 
tices are anisotropic. With NN forces only, they exhibit 
a vanishing modulus for compression applied in particu- 



lar planes rather than isotropically, but we expect them 
to exhibit many of the properties the two-dimensional 
lattices exhibit. Finally, we note that Maxwell's ideas 
can be applied to spin systems such as the Heisenberg 
anti-ferromagnet on the kagome lattice [4!| , and the 
possibility of unusual edge states in them is intriguing. 
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I 

Appendix A: Derivation of equilibrium, compatibility, and dynamical matrices. 



In this section, we will provide some of the details for calculating the equilibrium and compatibility matrices H and 
C and the Dynamical matrix D under periodic boundary conditions. The equilibrium matrix H relates the Nb = dN 
dimensional vector of bond tensions t to the dN dimensional vector f of forces at the nodes via H • t = f . We label 
unit cells by the their position vectors 1 = l\e\ + £2^2, where l\ and I2 are integers and 

2ix{n — I) 

e„ = (cos0„,sm0„); 0„ = (AI) 

are the primitive lattice vectors. (Here we take the lattice spacing to be one, i.e., 2a cos a I.) Each unit cell has 
three sites, labeled I, 2, and 3 as shown in Fig. [7j so we label the 3N C sites by (1, a), where a = I, 2, 3 and N c = N/3 
is the number of unit cells. Forces at the nodes are denoted by fi )0 and their 6N C = 2N components as f\. a .i, where 
i = x,y. They have to be balanced by stretching forces in the springs located on the bonds. There are six bonds per 
cell, which we can take to be the six bonds in the distorted hexagon shown in Fig. oriented parallel to the unit 
vectors 

( Tl — f )7T 

b„(a) = (cosV^sinV-n); ip n = - — + (-l)"a (A2) 

for n = 1, 6, where a is the twist angle of the twisted kagome lattice. When a = 0, these vector reduce to the edge 
vectors of a symmetric hexagon with b p (0) = — b p+3 (0). 

The 6N C = 2N bonds, labeled with 1, n, are occupied by springs under tension t\^ n that exert pulling (for t\ >n > 0) 
forces, direct along the bond vectors ±b„(o;), on nodes. We define vectors ti >n = ii jn b„(a;). Because we have periodic 
boundary conditions, we can express both fj „ and ti„ in terms of their Fourier components, 

(A3) 

ti,„ = -^=> V q l t q ,„. (A4) 



The equilibrium matrix is diagonal in q, so it breaks up into independent 6x6 blocks for each of the N c = N x N y 
independent vectors q. For simplicity, we consider only the case N x — N y for which N c = ATj and for which the 
number of sites is N = 37V,?. 




1. The Equilibrium Matrix 

There are four bonds incident on each site. The equations relating fi iU to ti jra follow from Fig. [7] 

fl,l = U,l — tl,6 + tl+e 3 ,4 — tl+e 3 ,3 (A5) 
fl,2 = tl,3 — ti2 + tl+ei,6 — tl+ei,5 (A6) 
fl,3 = tl,5 — t^ 4 + ti +e2 ,2 — tl + e 2 ,l- (A7) 
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FIG. 7: Watch this space 



We map (a, i) to a single index m with (a, x) — > 2a — 1 and (a, y) — > 2a and then Fourier transform. The result is 
that H breaks up into independent 6x6 blocks H nm (q) for each q. This leads to 



where 



H(q) = 



/ b hx (a) 
t>i, y (a) 
-&2,*(a) 
-6 2 ,2/(a) 

-e iq ' e2 6 2 ,*(a) e**«*&i, x (a) 

\-e^b 2 , y {a) e^h.yia) 



/«(q) = fln,m(q)tm(q)) 



-e iq ' e3 6 3 , K (a) e lq - e3 6 4ia; (a) 

-e lqe3 & 3 !y(a) e 4qe3 6 4 ^(a) 
&3,*(a) 
h, y (a) 

&5,z(a) 



(A8) 



"6,x 



(a) \ 



-e^b 5 , x (a) e^b 6 , x {a) 
-e iqei 6 5j2/ (a) e^beM 
-b iiX (a) 







bb, v {a) 



-bi, v {a) 







(A9) 



The null space of this matrix is easily calculated with the aid of Mathematica. 

Case 1, a = 0: This is the untwisted kagome lattice. When q ^ 0, the null space of H(q) is empty unless q is 
perpendicular to one of the primitive lattice vectors: 



q = q±, n e±,n = q(— sin <p n , cos n = 1,2,3, 

and q • e„ = qu n = 0, in which case there is a single vector, 



*i(q) = VNcS m ,,0(1,0,0, e"^- 163 ' 



ta(q) = V^c^ M , 2 ,o(0,l,0,0,( 



— iq±,2 ei-e. 



*s(q) 



iV c 5 ?||3i0 (0,0,l,0,0,e-^ 



,0,0) 

• 2 ,o) 



,3 e2-ej_ j3 



(A10) 

(All) 
(A12) 
(A13) 



where N c is the number of cells, for each n and q± : n in the null space. When q = 0, there are three vectors, 
(1,0,0,1,0,0), (0,1,0,0,1,0), and (0,0,1,0,0,1) in the null space. These are the q = limits of h(q), t 2 (q), and 
t 3 (q), respectively. We can, therefore, take the set of null-space vectors of H to be ti(q), i 2 (q), and t 2 (q) for the N x 
values of q, which include q — 0. Linear combinations of these vectors are also null-space vectors, and we can use 
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them to construct vectors confined to a single line of bonds. For example the vector 



*i(i) = j^E^Vq) = to,,o,o,o,* Isi _ V 3 /2) o,o) 



(A14) 



corresponds to a state in which springs on all 1-bond in cells with centers at 1 = ^ei for any l\ and on all 4-bonds in 
cells with centers at 1 = l\B\ +e 3 = {l\ — (1/2), — V3/2) are all under the same tension. This is a state of self-stress in 
which all bonds along a particular straight line parallel to the ex-axis are under tension. Similar states of self-stress 
for any of the 3N X lines of bonds can be constructed. 

Case 2, a > 0: This is the twisted lattice. The null space of H(q) is empty for all q > 0, and it contains only two 
vectors when q = 0. These correspond to two states of self-stress in which the stress the stress on the six bonds take 
on both positive and negative values that are identical in all unit cells. 



2. The Compatibility Matrix 



The compatibility matrix C relates the dN displacements d to the iVf, stretches si, a . The displacements are labeled 
the same way as the forces, i.e. as di j0 and the stretches in the same way as the bond tensions, i.e. as si ; „. The 
equations relating d\ a and si ; „ = — si, n b n (a) are 



dl-e 2 


3 - 


di,i 


= Sl,i = 


-si,ibi(a 


dl,2 


di 


-e 2 ,3 


= s l,2 = 


-si,2b2(a 


dl_e 3 


l — 


dl,2 


= s l,3 = 


-si,3b3(a 


dl,3 - 


di 


-e 3 ,l 


= Sl i4 = 


— si,4b4(a 


di_ ei 


2 — 


dl,3 


= si i5 = 


-«i,5b5(a 


di,i- 


di 


-ei,2 


= si i6 = 


-si i6 b 6 (a 



(A15) 



Fourier transforming and using the same (n, m) notation as for the equilibrium matrix, we obtain 



C mn (q)d„(q) = s m (q), 



(A16) 



where C mn (q) = H nm (— q) = Hj ran (q) as required by the constraint C = H 2 
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3. The Dynamical Matrix 

Finally, the dynamical matrix is D(q) = fcH(q) • C(q), whose components are 

1, 



D 12 

D 13 

D u 
D 15 
D 16 
D22 

D 23 

D 2i 

D 2 ,5 
D 2 ,6 

D 33 
D 34 

D 3fi 
D44 



£4,5 
D 5 , 5 

D 5 ,6 

D 6 ,e 



q) = -[cos(2a)+4] 

q) 
q) 



q) 

q) 
q) 

q) 

q) 

q) 

q) 
q) 

q) 

q) 
q) 
q) 

q) 



= — i\/3cos(2a) 



DIM) 



q) 

q) = D 5 * 1 (q) = -2 e -^, C( ^ (n)( , )S (^_ 

q) 
q) 



sin 2 Q - aj - e sin 2 (a + ^) 
(g-2a) + -c ^ ^cos(2a + -) 



DJ^q) = sin(2a) sin (|-) 

2 sin 2 (a) + cos 2 ^— — a) + cos 2 (a + — ) 

(*-*■) 



q) = ^32 (q) 



1 -<(*+ 
2 e 



^ 2 (q) 



cos' 



e i9a; cos 2 



£>5 i2 (q) = ie - '^ 9 " sin(2a) sin 
^ 2 (q) = -2 e -^^sin 2 ( a )cos(|) 
2 sin 2 Q - a) + 2 sin 2 (a + ^) 

£>43(q) = o 

■i^-*) . 2/7T N 



As, 3 (q) 



V 



e i9x sin 2 



^6, 3 (q) = -2 e ^ C0S (J- 2a )"2 

2 cos 2 Q - a) + 2 cos 2 (a + ^) 

^ 4 (q) = -^ cos (J- 2a) -\e 

DiM = -e^ 2 2 ^cos 2 ( 
sin 2 (J^ - aj + sin 2 (a + ^) + 2 cos 2 (a) 

^6, 5 (q) = 5 cos - 2a) + 1 cos (2a + |) 
2 sin 2 (a) + cos 2 ^ — a) + cos 2 + ^) 



cos 



cos 



a ] - e tq * cos 



( 2a+ 6 



(A17) 



Its eigenvalue spectrum can easily be calculated with the aid of Mathematica. The results are shown in Fig. 3 of the 
main text. 



Appendix B: Surface Modes 



Surface modes decay exponentially into the bulk. They are characterized by a wavevector q|| parallel to the plane 
of the surface, their frequency w s (q|i), their decay length /(q||) perpendicular to the surface. In lattices with nearest- 
neighbor forces only, w s (q||) and Z(q||) can be determined by setting q = q|| + q_i_, where q is the wavevector that 
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appears in the dynamical matrix and qj^ = q±e± is the component of q perpendicular to the surface (ej_ is the unit 
vector perpendicular to the surface), setting q± = iZ , requiring 

det[u; 2 /- LHqn^r 1 )] = 0, (Bl) 

where I is the unit matrix, and that the equation of motion of the surface layer be satisfied. In the case of surface 
modes of with zero frequency, the relation det[D(q|| , iZ -1 )] = determines Z as a function of q||. The surface modes 
of the twisted kagome lattice have zero frequency, so we need only solve this equation. We consider only the case in 
which the surface is perpendicular to the y-direction and set q x = q and r = e l ^ qy l 2 = e~^ , K 2l \ We find 

det[D(«, r)] = ^a(q) [r + r~ x - 2g(q)] [r + r" 1 - 2g*(q)] , (B2) 



where 



a(q) = 2 [2 - cosq + (1 - 2 cos g) cos(4a)] (B3) 

cos(g/2)[3 — 2 cosp + cospcos(4a)] + 2iv3sin 3 (3(7/2) sin(4a) 
9 ^ q ' ~ 2 - cos q + (1 - 2 cos q) cos(4a) 



Solving for r in Eq. (|B2j) and to a lattice spacing of 2a cos a, we obtain 

1 



k y = I 1 = — -= log[— g(2aqcosa) + y / g(2aqcosa) 2 — 1]. (B5) 

av3cosa 

This is the function that is plotted in Fig. 5(b) in the text. Note that k y has both a real and an imaginary part, but 
that in the limit q — > 0, k y = q. 
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